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ABSTRACT 



Aims. The non-conservative gravitational potential of barred galaxies, or of any other non-axisymmetric structure, produces a loss of 
energy in infalling particles of the intergalactic medium into the galaxy, which are trapped in its potential. This dynamical friction can 
contribute towards increasing the total mass of barred galaxies. 

Methods. Analytical calculations of the energy loss are carried out using the orbits of the particles derived numerically. Theoretical 
predictions are compared with observations through the statistical analysis of the rotation curves of barred and non-barred galaxies, 
either in cluster or field galaxies. 

Results. There is a net effect of accretion, but it is normally very low in relative terms. It is only significant (> 10% of the total mass 
of the dark matter halo in the life of the galaxy) if the density of the intergalactic medium is higher than ~ 3 x 10 13 M /Mpc 3 (or 
considerably lower in cases of motions of the galaxies close to the IGM average motion, or perpendicular to the plane of the galaxy, 
or when the halo mass is low). 

Data on rotation curves do not show clear trends towards higher halo mass for barred galaxies, only slight trends for early-type spiral 
galaxies. In any case, the statistical uncertainties are limited to the detection of differences in masses >a;20%, so the effect of bar- 
driven injection of intergalactic matter into galactic halos might be present with relative contributions to the average mass of these 
barred galaxies lower than 20%. 
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1. Introduction 

Intergalactic matter (IGM) is still a field of astrophysics with 
many open questions and few answers. There are either theoreti- 
cal arguments or interpretations of observations that suggest the 
existence of an IGM, although the definitive details of impor- 
tant parameters such as the mean density, whether continuous 
.J^i | or in the form of clouds, etc. are not clear yet. Cold dark mat- 
, ter models predict substructure within galactic- and cluster-mass 
' halos that form in a hierarchical accretion Universe (Klypin et 
?3 \ al. 1999; Moore et al. 1999) that could constitute some of the 
clumpy structures in the IGM. Indeed, stellar and cold gas in 
galaxies contribute 8(^)% of the total amount of the Big Bang- 
produced baryonic matter predictions (Bell et al. 2003); this 
implies a low overall efficiency of galaxy formation, and that 
the rest of the baryonic material must be some way away from 
the visible galaxies, possibly in the IGM. Some interstellar gas 
is ejected into the IGM by tidal interactions between galaxies 
(Morris & van den Bergh 1994), SN explosions, etc. Of course, 
the existence of intracluster gas in rich and irregular clusters, 
which have diffuse emission and X-ray emission associated with 
hot gas, is well established. Possibly, some high-velocity clouds 
(HVCs) observed in the radio are located outside our Galaxy, i.e. 
belonging to the IGM (Blitz et al. 1999; Braun & Burton 1999). 
There is also a large number of lines visible in optical at high 
redshifts or UV in low redshifts — the Lyman-ff forest — which 
are presumably caused by clouds along the line of sight of QSOs 
(Rauchl998). 

There is a need for a substantial accretion of low angular 
momentum material from the IGM into the galaxies (Fraternali 
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et al. 2007). Chemical evolution in the local Galactic disc (the 
G-dwarf problem, 9 Be abundance and others; Rocha-Pinto & 
Maciel 1996; Casuso & Beckman 1997) needs a continual but 
episodic infall of metal-poor gas that mixes slowly with the rest 
of the interstellar medium. Lopez-Corredoira et al. (1999) de- 
duce that the IGM mass must currently represent at least around 
one half of the total mass of the Local Group, given that the 
accretion rate, as inferred from chemical evolution, has not de- 
creased significantly during the disc lifetime. 

The existence of matter in the IGM should have some effects 
on the formation and evolution of the galaxies. For instance, it 
is a possible explanation for the formation of warps in spiral 
galaxies (Lopez-Corredoira et al. 2002a). It is also, as said, a 
likely explanation for the chemical evolution of galaxies (Casuso 
& Beckman 1997). Here, we study another effect: the accretion 
of this IGM matter due to dynamical friction. It is clear that the 
disc mass might be fed by the collisional friction of baryonic 
matter with disc gas. The enrichment of the halo by means of a 
gravitational interaction is instead the topic to be treated here. 
In particular, the component of the galaxy responsible of this 
gravitational interaction will be the bar. 

Bar-halo interactions have been studied for several reasons. 
For instance, as a mechanism of angular momentum transfer. 
It is well established that strong bars rotating in dense halos 
generally slow down as they lose angular momentum to the 
halo through dynamical friction (Debattista & Sellwood 1998, 
2000; Athanassoula 2003, 2005). This friction can be avoided 
or found to be anomalously weak in some circumstances and 
bar slowdown can be delayed for a period in a metastable state 
(Sellwood & Debattista 2006), although Sellwood & Debattista 
(2006) demonstrate that mild external, or internal, perturbations 
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quickly restore the usual frictional drag; it is therefore unlikely 
that a strong bar in a galaxy having a dense halo could rotate 
for a long period without friction. A more long-lasting effect to 
be considered is that the velocity dispersion of the halo parti- 
cles is high enough to stop the resonances from absorbing con- 
siderable amounts of angular momentum (Athanassoula 2003), 
hence prohibiting bar slowdown. Owing to this angular momen- 
tum transfer and self-consistent re-equilibration, strong realistic 
bars will modify the cusp profile, lowering the central densities 
within about 30% of the bar radius in a few bar orbits (Weinberg 
& Katz 2006). 

In this paper, another effect of the bar will be explored, also 
using the dynamical friction of this kind of gravitational interac- 
tion: the influence of the bar to produce the accretion of IGM 
matter (initially non-gravitationally linked to the galaxy) into 
the halo. That is, the interaction with the bar will produce the 
loss of energy of some particles that cross the galaxy and that 
will be trapped in it, since they do not have escape velocity, to 
form part of the galactic halo. This mechanism applies to both 
baryonic and non-baryonic matter because the friction is purely 
gravitational. In §|2] this loss of energy and the amount of ac- 
creted matter depending on IGM parameters for a Milky-Way- 
like galaxy is calculated semianalytically. In §|3] it is explained 
that other factors of the non-conservative potential, for example 
approaches to individual stars, produce a much lower effect that 
is totally negligible. In §|4] statistics with observational data are 
performed in order to see whether barred galaxies have larger 
masses. 



is conserved. 

Let us consider a monodimensional bar. This is an approx- 
imate description of real 3D-bars because the major axis is 
considerably larger than the other axis — otherwise, they would 
not be bars — as can be observed, for instance, in the galaxies 
NGC 1300, NGC 7479, NGC 4123, NGC 1433, NGC 4999, 
etc. (Sandage & Bedke 1994) or in our own Galaxy (Lopez- 
Corredoira et al. 2001, 2007). The mass density distribution will 
be A(l) as a function of the distance / to its centre. The azimuthal 
component of the gravitational force produced by this bar of ra- 
dius R{) on a particle of mass m will be: 



FJR,(f>',z) = -Gm 



J_ Ro U \ Z 2 + R2 + 



A{l)l sin 4>' 



(6) 



■l 2 -2Rlcos 0') 3/2 
If we assume A(l) = constant = A, the integral can be solved: 



FJ r '(R,<f>',z)] = 



-GmA sin <p' 
z 2 +R 2 (l-cos 2 4>') 



(7) 



Ri R {) R cos 4>' 
r+ + 

r+ r+ 



Rl RoRco&f^ 
r_ H i 



r + = \r' - Ho|; r_ = \r' + R \,R = R L 

Therefore, if we want to calculate the gain/loss of energy of 
the particle, we must just know its orbit r(t) and use expressions 
© and © to calculate AE. 



2. Variations of energy in a barred potential 

We analyse the variations of the energy in particles that follow 
orbits in certain kinds of non-conservative potentials: those in 
which the gravitational force varies with time. This happens, 
for instance, in a barred galaxy, because the rotation of a non- 
axisymmetric structure produces variations in the potential due 
to the bar with respect to an inertial frame. However, the rotating 
disc has a stationary potential since the axial distribution of mass 
does not change. If the bar rotates with angular velocity to, the 
variation in the energy of a particle along its path (orbit) is: 



A£= ( F(R,cf>- ajt,z)dr 

J orbit 



(i) 



; 

J orbit 



F(R, <p - Lot, z)(dReR + Rdfaf + dzk), 



in cylindrical coordinates, where F is the gravitational force pro- 
duced by the bar. If we change to the variable 4>' — cp - cot (that 
is, we adopt the frame that is moving with the bar): 

AE = I F(R, <f>',z)(dRe R + Rd<f>'e r + dzk + Rtodte#). (2) 

J orbit 

If the orbits of particles begin at infinite distance and finish 
at infinite distance: 



I 



F{r')dr' - (because it does not depend on time), (3) 



r=oo, orbit 



AE 



■£ 



FAR,<p',z)Rdt = uAJ z . 



(4) 



vbit 



The Jacobi integral 
H = E-a)J, 



2.1. Numerical calculation of orbits 

For our problem, we use semianalytical calculations. The solu- 
tion of eqs. 0}, $6fy is calculated numerically, and the orbit of 
each particle itself is not solved numerically as an outcome of 
a numerical simulation in which all particles interact with each 
other but as the numerical integration of its motion equation. 
Indeed, we can consider the gravitational potential of the in- 
falling particles as due only to the components of the galaxies 
and neglect the interaction due to other IGM particles. 

Analytically, we cannot calculate the orbit of a particle in a 
galactic potential; the integrals are too complex to find the solu- 
tion. But we can calculate r{t) simply by numerical calculation 
of the differential equation f = a ga i axv , the total acceleration pro- 
duced by the total potential of the galaxy. It is then straightfor- 
ward to obtain r'(f) (the position in the frame that is moving with 
the bar) from r(t) by using <p' - (p-cot. In this case, of course, for 
the calculation of the orbit, all the potentials are to be taken into 
account not only the potential of the bar. We can approximate 
this potential as the sum of several components: 



% galaxy ~ ^halo >" ^disc ^ba: 



(8) 



We give analytical expressions for these acceleration terms as 
integrals. All these integrals will be performed numerically with 
the Simpson algorithm, and the solution to the differential equa- 
tion f = ag a iaxy will be carried out with a variable step depending 
on the acceleration: 



At = 



50s 



a(m/s 2 ) + 10 



•12 ' 



(9) 



The Jacobi integral H is conserved for each orbit within vari- 
ations of ~ 0.1% . The accuracy in the orbit calculation was also 
checked by testing that the orbit is approximately the same when 



(5) decreasing this step At by a factor of 50 (At 



Is 



o(m/s 2 )+10- 



r ): the 
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difference in the position of both calculations after a travelled 
path of more than 1000 kpc (once the particle fell from a dis- 
tance ro towards the centre and, after that, it crosses the shell at 
r = ro; with ro = 400 - 500 kpc) is ~ 0.4 kpc on average. These 
errors are not important for the present calculations because we 
do not require too much accuracy in the orbit; it is only used 
for the path integral of eq. (@). Moreover, a small random error 
added to the orbit is equivalent to a small random variation of the 
initial conditions, which in comparison with the range of impact 
parameters (0-200 kpc) is negligible. 

2.1.1. Halo potential 

We approximate the halo as spherical, i.e. we neglect the ellip- 
ticity: 



ahaioW 



[4nG £ dt t 2 p hah (t)]i 




(10) 



We then use a truncated flat model by Battaglia et al. (2005): 
M halo a 2 



Phalo(r) = 



4n r 2 (r 2 + a 2 ) 3/2 ' 



(11) 



2.1.2. Disc potential 

For the disc, we use a typical exponential function considering 
its thickness to be negligible. In order to reduce the calculation to 
only one integral, which reduces the calculation time very con- 
siderably, we express it in elliptic integrals [Binney & Tremaine 
1987, eq. (2-142b)]: 



a d i S c[r(R,(f>,z)] = — 7rrr e R ^ — «> 



dR 

■2n 



dz 



dTT I 

o Jo 



da- 



G<t(T) 



Fig. 1. Modulus of the acceleration in a monodimensional bar 
potential with the parameters of {32.1.41 Logarithmic scale. 
Range of x and y between -20 kpc and 20 kpc. 



where E(k) are the complete elliptic integrals of the second kind, 
E(k) 



yjz 2 +R 2 + T 2 - 2RT cos(0 - a) 



(12) 
(13) 



- z^L 

where 
KR, T,z) 



dT K(k)ko-(T) VF, 



ART 



]z 2 + (R + T) 2 ' 
and K(k) are the complete elliptic integrals of the first kind 
dt 



(14) 



K(k) 



-f 

Jo 



V(i-? 2 )(i-£ 2 f 2 ) 

The two derivatives of d> are: 



<90disc 



f 

Jo 



X 



OR R 3 ' 2 jo 
It k 2 



dT 



(15) 



(16) 



K(k) 



R 



-- + — \E(k) 
A\\-k 2 )\R T RT ) 



kcr{T) Vf 



(Binney & Tremaine 1987, eq. [2-146]), and 
Gz f°° ,„k 3 E(k)o-(T) 



"disc 



dz 



2R 3 ' 2 



r°° k 3 E 



k 2 Vr 



(17) 



i - k 2 t 2 



J0 ' i-f 
We assume an exponential disc such that: 

R 



cr(K) = cr exp|- — 



2.1.3. Bar potential 



(18) 



(19) 



For the bar, as said, we adopt a monodimensional structure ro- 
tating with angular velocity to; 



«bar = -G I dl 

J-Ro 



A(l)[r - l[cos(a>t)i + sm(cot)j]] 
(z 2 +R 2 + l 2 - 2Rlcos((p - ojt)) 3 ! 2 ' 



(20) 



Other bar potentials could be used (e.g., Athanassoula 1992) but 
this one is used for its simplicity. The IGM particles will never 
cross this ID structure, because there is zero probability that the 
path by chance crosses the line of the bar, so this potential will 
not give problems of divergence. The potential also behaves well 
at large distance, as shown in Fig. [T] the acceleration falling as 



2.1 .4. Parameters for our numerical calculations 

The parameters that we used are those more or less typical for 
Milky Way-like galaxies. 

A total mass for the halo M ha i = 1 -2 x 10 12 M Q , and a = 105 
kpc (Battaglia et al. 2005). The mass within a radius r would be 
M(r) = Mhaio ?4 f ■ 96.7% of the halo mass is within r < 400 
kpc, so it can be approximated for r > 400 kpc as a Keplerian 
potential of mass equal to Mhaio (plus the mass of the disc and 
the bar). 

For the disc, we adopt ctq = 1 kg/m 2 (in agreement with 
the solar neighbourhood measurements by Kuijken & Gilmore 
1989) and Iir = 3.4 kpc (a typical value in a Milky Way-like 
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galaxy; Lopez-Corredoira et al. 2002b). The corresponding total 
mass ofthediscis3.5xl0 10 M o . 

For the bar, we assume a total mass of 1.7 x 10 10 M 
(Sevenster et al. 1999), more or less in agreement with Kuijken 
& Dubinski's (1995) self-consistent models of the Galaxy. This 
gives a constant distribution A{1) — 1.4 x 10 20 kg/m for a radius 
Ro = 4 kpc (estimated value of the long-bar radius in our Galaxy, 
Lopez-Corredoira et al. 2001, 2007). This mass of the long bar 
is indeed not well known (Lopez-Corredoira et al. 2007) and 
the bar of Sevenster et al. (1999) is 35% shorter (the axial ratio 
difference is not so important, as we will see below) than Lopez- 
Corredoira et al. (2007) bar, so we are not setting the parameters 
of the bar accurately, but as a more or less correct order of mag- 
nitude (see further discussion in {32.6l of what would happen with 
more or less massive bars). 

The angular rotation velocity, to, is taken as 1.9 x 10~ 15 s _1 
(as measured in the Milky Way; Debattista et al. 2002). Most 
strong bars in real galaxies appear to rotate rapidly (Aguerri et 
al. 2003), so it is not expected that this velocity will slow down 
too much. Although the dynamical halo-bar friction predicts a 
quick slowdown in this velocity (Debattista & Sellwood 1998, 
2000; Athanassoula 2003, 2005), there might be a metastable 
state in which the bar-halo also remains a long time without 
friction (Sellwood & Debattista 2006), or a state in which the 
velocity dispersion of the halo particles is high, so as to stop 
the resonances from absorbing considerable amounts of angu- 
lar momentum (Athanassoula 2003), thus prohibiting bar slow 
down. Whatever the explanation, the fact of fast bars in most 
barred galaxies remains, and we therefore do not consider any 
friction or change in co with time. 

With these parameters, the stationary (Lagrangian) points of 
the Galactic potential are at: L1,L2: <p' L — 0, n, z — 0, Ri = 5.25 
kpc; L4,L5: <f>' L = tt/2, 3tt/2, z = 0, R L = 4.59 kpc (for more in- 
formation about the Lagrangian points in a barred potential see 
Pfenninger 1990 or Athanassoula 1992), whose Jacobi integrals 
are: H = -1.41 x 10 11 J and -1.29 x 10 11 J respectively. The fact 
that the Lagrangian points L4,L5 are 13% closer to the centre 
than the Lagrangian points L1,L2 is expected for average strong 
bars; the difference is typically 6% for average bar models ac- 
cording to Athanassoula (1992) but it may be as large as 20% 
for the strongest bars. The strength of the bar can be compared 
with another parameter defined by Athanassoula (1992): half of 
the quadrupole moment along its major axis 

Qm = \M bM Rl (21) 

This expression is derived from eq. (5) of Athanassoula with 
a - Rq, b — (infinite axial ratio a/b), n — (constant density 
along the bar). With our parameters we get Q m - 5.4 x 10 4 (in 
units 1O 6 M kpc 2 , as in Athanassoula 1992). This is again an in- 
dication that the strength of our bar is a typical one within strong 
bars (weak bars have Q m ~ 10 4 , very strong bars have Q m ~ 10 5 ; 
Athanassoula 1992). It is also remarkable that the infinite ax- 
ial ratio a/b in a monodimensional bar does not significantly 
alter the effects of a bar with finite but high axial ratios since 
Q,„ cc (a 2 - b 2 ) (Athanassoula 1992, eq. 5) and (a 2 - b 2 ) » a 2 for 
high a/b, either finite or infinite. 

2.1 .5. Initial conditions of the IGM particles 

We assume now that our galaxy is moving with respect to the 
IGM (or the IGM is moving with respect to the galaxy). This 
hypothesis was already used in Lopez-Corredoira et al. (2002a) 
as an explanation to produce warps in disc galaxies and now will 



be explored just to calculate the fraction of the particles that are 
gravitationally trapped in the galactic potential due to a negative 
energy variation. Basically, the model consists of a constant den- 
sity (pea) medium which has at infinite distance from the galaxy 
an average velocity v,*, with respect to the centre of galaxy, and 
an angle of with respect to the plane of the galaxy equal to 6U . 
We neglect the effect of the dispersion of velocities with respect 
to this average velocity. 

We randomly distribute the values of b (impact parameter 
with respect to the centre of the galaxy), <pb (azimuthal angle 
in cylindrical coordinates with the vector v«, as vertical axis) 
following the conditions for a homogeneous distribution, equiv- 
alent to do integrals in the way J d(pb J """ db b (...). We set 
bmax = 200 kpc. That is, we neglect the contribution of particles 
whose impact parameter is larger than 200 kpc; this is justified 
because, as we shall see in following subsections, the particles 
which lose more energy have lower value of b, except in cases 
with low Voo. 

Since an infinite distance cannot be introduced in the com- 
puter (integration of the orbit would take an infinite time), we 
must begin the integration from a finite distance, which we 
choose to be at 400 kpc< ro < 500 kpc. We distribute our parti- 
cles randomly at a shell with this range of radii. The distribution 
with different ro is useful to integrate over different positions of 
the bar when the particle reaches the galactic disc; that is, it is 
equivalent to doing an integration with different initial bar an- 
gles «(fo). The time to cross 100 kpc at around 70 km/s (the 
velocity of the particle at radius ro assuming v M = 50 kpm/s) is 
1 .4 Gyr, so doing the integration between 400 kpc< ro < 500 
kpc is equivalent to doing an integration over a time of 1 .4 Gyr, 
equivalent to 14 cycles of the bar (one cycle/100 Myr); more- 
over, there is an added dispersion in the position of the bar when 
the particle reaches the centre of the galaxy because the different 
values of </>/, and b get different infall times to reach the centre, 
so in practice we are integrating over a range larger than 1 .4 Gyr. 
This ensures that all angles of the bar are covered more or less 
uniformly (departures from uniformity are negligible and less 
important than Poissonian or other fluctuations). In this way, we 
are doing a Monte Carlo simulation over the variables b, 4>b and 
ro- 

At the distance of 400-500 kpc the non-monopolar terms of 
the gravitational potential are negligible and almost the entire 
mass of the galaxy is within this radius, so the potential can be 
approximated as one produced by a point-like source in the cen- 
tre of the galaxy with its total mass. Nonetheless, this potential 
is not totally negligible at these distances, so we must correct 
the initial conditions with respect to the conditions at infinite 
distance: 



r o = - Jro - r 2 mp e Voa + r imp e bm , (22) 

where e Va> is the unit vector parallel to the initial velocity at in- 
finity, and eb m is the vector perpendicular to e Voo in the plane per- 
pendicular to the one defined by ro and v^,, i.e. the direction of 
the impact parameter, pointing towards an increase in b. 

vo = vo cos(or - ao)e Vx - vq sin(a - ao)e b„ ■ (23) 

The parameters Vo, r- imp and a, ao are derived from the analyt- 
ical expressions of a hyperbolic orbit (see Lopez-Corredoira et 
al. 2002a; subsect. 3.3) and taking into account that the angular 
momentum is constant along the orbit: 

/ 2GM(r ) ; 2 

vo = \ + vto, (24) 
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n mp = r \ sin^S- cos (-l/e)]|, 



(25) 



(3 = cos 



v. =50 Km/s, 8 =45 , M =1.2 x 10 M. , 1x200 kpc 



GM(r ) e 2 - 1 ' 



1 + 



GM(r )) 



(26) 



(27) 



(28) 



ao = — I- cos 1 1 — - I (defined to be between n/2 and n), (29) 
2 \ r 



a = 7r - sin 



(denned to be between 7r/2 and n). (30) 



The results of the orbit with these analytical expressions 
were compared with the results of numerical calculations, and 
they agree. By applying these formulae, we do not need to do 
numerical calculations beyond 400-500 kpc. 

2.2. Histogram of energy variations 

The only free parameters for the calculation of the energy varia- 
tion of each particle once the galactic parameters are established 
are p M , f«, and 9^. We could also consider Mhaio as a third free 
parameter instead of a fixed value of Mh a i = 1.2 X 10 12 M , 
since the total mass of the halo might be variable with the time, 
precisely because of the accreting mechanism proposed in this 
paper: the loss of energy of some particles which are trapped 
when this loss of energy is higher than 



IT 1 2 



(31) 



As an example, let us set v x = 50 km/s (a typical value for 
the peculiar motion of a galaxy) and 6 X = 45°, with Mh a i 
1 .2 x 10 12 M . The results for 10000 particles are in Figs, [3] 
and |4] which show that the loss of energy is high enough to be 
trapped only for b < 100 kpc or -2.2 x 10 11 < H < 1.2 x 10 11 
J [for comparison, the Jacobi integrals of the Lagrangian points 
are H = -1.41 x 10 11 J and -1.29 x 10 11 J], and only 426 of 
the 10000 particles lose an energy higher than E^. Equations 
© and (0 were applied to calculated AE for each particle once 
it moves away from the galaxy. There are, as observed in Fig. 
[2] practically the same number of particles which gain E dE as 
those which lose E dE. 

Let 7?2oo be the ratio of particles within b < 200 kpc which 
are trapped by the halo. According to the above, this is /?200 = 
= 0.0426 + 0.0021. Hence, the accretion ratio will be: 



dM t 



halo 



dt 



f d4> b f 
Jo Jo 



■200 kpc 



dbb 



(32) 



1.9 x 10 2 Voo(kg/m 3 )voo(km/s)i?200 M /yr 



= 4.0(±0.2) x 10 24 P oo(kg/m 3 ) M /yr 




Fig. 2. Histogram of distribution of energy gain/loss per unit 
mass (kg). 
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Fig. 3. Energy gain/loss in the simulation of particles with dif- 
ferent impact parameter (b). 



2.2.1. Dependence on the initial infall angle 

The dependence on the angle is illustrated in Fig. [5] The R200 
ratios are: 0.0441 (+0.0021) for e m = 5°, 0.0436(±0.0021) for 
Oct, = 25°, 0.0536(+0.0023) for ft*, = 65°, and 0.1563(+0.0040) 
for 9oo = 85°. R20Q does not change appreciably except for 9^ = 
85°, which is 3-4- times higher. The polar accretion seems to 
favour the mechanism. The dependence is roughly described by 
a law: 



RiooiOoo; voo = 50 Km/s, M Mo = 1.2 x 10 12 M ) 
« 0.043 + 0.20 exp[-0.12[90-6Udeg.)]] 



(33) 



Since we are interested in calculating the order of magnitude, 
we can disregard this dependence and leave for our remaining 
calculations a value of 9^ =45°. 
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Fig. 4. Energy gain/loss in the simulation of particles as a func- 
tion of its Jacobi integral. For comparison, the Jacobi integrals 
for the Lagrangian points are H = -1.41x10" J and -1.29xl0 10 
J. 



v w =50 Km/s, M Mo =1.2 x lo' 2 M sun , b<200 kpc, AE<0 




Fig. 5. Histogram of distribution of energy loss as a function of 
the initial infall angle. 



2.2.2. Dependence on the initial velocity 

However, the dependence on initial velocity is high. The his- 
tograms in Fig.[6]show it. The ratios /?200 are (note that E x varies 
withvoo): > 0.4966 forvoo = 1 Km/s, 0.25 11 (±0.0050) for v M = 
30 Km/s, 0.01 18(±0.001 1) for v x = 70 Km/s, 0.0035(±0.0006) 
for v M = 90 Km/s, and 0.0012(±0.0003) for = 110 km/s. 
The ratio of higher energy loss is maximum for lower initial 
velocities, both in terms of absolute energy and relative energy 
with respect to E^. The explanation for this is multiple. First, for 
higher velocities the particles go very fast when they reach the 
bar, and since AE in equation (0]i is proportional to the time the 
particle spends near the bar, the loss of energy is lower. Second, 
the capacity to attract towards the centre particles with high b is 
larger for lower values of the initial velocity. Moreover, E x is 
higher for higher initial velocities, so the rate of energy loss is 



Fig. 6. Histogram of distribution of energy loss as a function of 
the initial velocity. 

lower too. For the cases Voo < 30 km/s, it is observed that some 
particles with b > 200 kpc also lose AE < therefore, the 

trapped mass is even higher than the ratio given in eq. d32b . The 
dependence is roughly described by the law: 

#20o(Vco 30 km/s; M = 45°, M halo = 1 .2 x 10 12 M ) (34) 

~ 0.82exp[-0.06vco(km/s)]. 



2.2.3. Dependence on the halo mass 

The effect of the variation of halo mass is similar to the ef- 
fect of the variation in the initial velocity because the amount 
of mass is proportional to the acceleration, and consequently 
is related to the velocity of the particle. Fig. [7] illustrates this. 
A maximum accretion is obtained for the largest mass. Lower 
masses have lower capacities to trap particles with high b. The 
values of R 20 o obtained are: 0.0049(±0.0007) for M halo = 0, 
0.0092(±0.0010) for M ha i D = 0.2 x 10 12 M , 0.0220(±0.0015) 
for M halo = 0.4 x 10 12 M , 0.0299(±0.0017) for M halo = 
0.6 x 10 12 M Q , 0.0362(±0.0019) for M ha , = 0.8 x 10 12 M , 
and 0.0396(±0.0020) for M ha , = 1.0 x 10 12 M . It is clear that 
it is proportional to the halo mass, roughly according to: 

#20o(M ha i ; 0„ = 45°, v., = 50 Km/s) (35) 

* 0.0062 + 0.034M ha , o (10 12 M ) 

2.3. Other ways of accretion and the escape of halo particles 

The gain or loss of mass in the halo due to other mechanisms will 
not be considered in this paper, but only the mechanism related 
to the dynamical friction with the bar. There are other mech- 
anisms of accretion, such as friction with the disc of particles 
or clouds that cross it. There are also ways for the halo to lose 
mass, such as galaxy stripping in interactions with other galaxies 
or the IGM, or the escape of particles when they reach a velocity 
larger than escape velocity due to multiple interactions. We ne- 
glect these other mechanisms because we want to consider only 
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2.4. Evolution of the halo mass 



v. iit .=50 Km/s, B. n( =45 , b<200 kpc, AE<0 




Fig. 7. Histogram of distribution of energy loss as a function of 
halo mass. 



Let us suppose that the dependence on the initial angle, initial 
velocity and the halo mass are separable. This is not strictly cor- 
rect, but may be appropriate as a rough estimation . In this case, 
a generalization of the expression d32b will be: 



—P^ = 2.1 X 10 2 V(kg/m 3 )v ra (km/s) x (37) 
at 



exp[-0.06vc»(km/s)][l + 4.7exp[-0.12[90 - 6Udeg.)]]] 
x[l + 5.5 x l(T I2 M halo (M )] M /yr. 

Here we are assuming that the disc and the bar do not change 
with time. The solution of this differential equation will be: 

MhaioC*) = -1.8xlO 1I M +[1.8xlO I1 M +M ha i o (f = 0)]xe £ ,(38) 

£=1.2x 10 12 poo(kg/m 3 )voo(km/s) x exp [-0.06voo(km/s)] 



the effect of the bar in increasing the halo mass independently 
of other effects. The reality is more complex, of course, but here 
we just consider a toy model for rough estimates of the net effect 
of an increase in the accretion of mass in barred galaxies with 
respect to non-barred galaxies. 

The non-conservative potential of the bar could also produce 
an increment of the energy of the halo particles and these could 
reach the escape velocity producing some kind of "evaporation" 
of the halo. The mechanism is the same as the one explained in 
previous subsections but the initial conditions will be different 
from those explained in 32.1.51 Indeed, the fact that we get some 
net accretion of some IGM particles arises because some of them 
go through the centre of the galaxy, but the virialized halo has 
totally different initial conditions. 

Assuming that the halo particles have an initial velocity cor- 
responding to a circular velocity plus isotropic dispersion veloc- 
ities with a component in the radial direction [given the isotropy, 
the dispersion in other directions will be the same; the isotropic 
representation, while more or less correct in the inner halo, is not 
correct in the outer halo (Dehnen & McLaughlin 2005), but we 
consider it here as a rough approximation] 



<7>(r) oc [p(r)r a ] 



.a-,1/3 



x[l +5.5 exp [-0.12[90-Mdeg.)]]] xt(yr), 

and the total accretion in the lifetime of the Galaxy (we take 
t — 1.4 x 10 10 yr) due to this mechanism is 

AM halo = [1.8 X 1O U M + M halo (f = 0)] x [e E - 1], (39) 
E = 1.7 x 10 22 poo(kg/m 3 )v oo (km/s) x exp [-0.06v M (km/s)] 
x[l + 5.5 exp [-0.12[90 - 0»(deg.)]]]. 

The necessary conditions for a fraction F of the halo mass to 
be accreted by this mechanism [Mh a i (f = 0) = (1 — F)Mhaio(0> 
AM = FMh a \ (t); assuming there are no other mechanisms of 
accretion] are: 



1 + 



1.8 x 10 U M P 



M ha i (f) 



x 1 + 



1 



(40) 



(36) 



Particularly, for M halo (f) = 1 .2 x 10 12 M , 6» M = 45° and v*, = 50 
km/s, 



a = 35/18 (Dehnen & McLaughlin 2005), in agreement with 
simulations and analytical solutions of self-gravitating dynami- 
cal equilibriums; with the halo density of eq. ( ITTb , and cr r (r = 
a) = v c {a)l V3 (Dehnen & McLaughlin 2005), where v c (r) is the 
circular velocity, we can perform again the numerical calcula- 
tion of the orbits and the increment of energy according to eqs. 
(0|l and (|6]l. The result with 10000 halo particles whose initial 
conditions are distributed randomly according to the probability 
densities of positions and velocities including the dispersions, 
during 14 Gyr orbit calculation, gives the result that none of 
these particles gets an increase in energy AE > c ^' o) , where ro 
is the initial distance from the centre of the galaxy. This means 
that the increase in the halo evaporation due to the bar-halo in- 
teraction is negligible: ^ <~ 10~ 4 during the whole life of the 
galaxy. 



\ exp[4.3xl0 22 pc O (kg/m 3 )]-l/ 

This expression has sense only for pc^kg/m 3 ) < 4.7 x 10~ 23 
kg/m 3 (F < 1), for higher densities the mass of the halo should 
be larger. In Fig. [8] the values of F for different densities are plot- 
ted. Densities of the IGM around 10~ 24 kg/m 3 (including bary- 
onic and non-baryonic matter) are able to produce a significant 
amount of accretion by means of this mechanism. 

In order to have at least 10% of the total halo matter accreted 
by this mechanism (F = 0.1) it is necessary that (we keep 6^ - 
45°,M halo (f) = 1.2xlO 12 M ) 

F>0.1o v oo (km/s)e-° 06, ''» (km/s) poo(kg/m 3 ) > 5.3 x 10~ 24 (42) 
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Fig. 8. Fraction of halo mass accreted by means of the mecha- 
nism proposed in this paper. 



2.5. Considerations about the range of possible values of the 
parameters 

2.5.1. Density of the IGM 

The expression d42l for F > 0.1 indicates that for 6 X = 45°, 
MhaloW = 1.2 X 10 12 M , Voo = 50 km/s we need a density of 
2 x 10~ 24 kg/m 3 = 3 x 10 13 M /Mpc 3 = 10 2 /i 2 p critical . This 
is a high density and may be reached by the intracluster gas 
in some groups or clusters of galaxies, especially in the rich- 
est ones, which can have a density much higher than this (in- 
cluding dark matter). The typical IGM densities in clusters are 
between 10~ 25 and 10~ 23 kg/m 3 (Roussel et al. 2000). However, 
the relative velocities in such environments are also higher, thus 
decreasing the desired effect (see the next subsection). For cases 
such as the Local Group this density is too high. The average 
IGM density of the Local Group is <~ 10~ 25 kg/m 3 (assuming 
a total IGM mass of <~ 2 x 10 12 M G [Lopez-Corredoira et al. 
1999] in a volume of ~ 1 Mpc 3 ). 

2.5.2. Velocity of the IGM 

The product Vooe~ 006v °° is maximum for — 17 km/s, at which a 
density 9 xl0~ 25 kg/m 3 produces according to eq. d42l a F = 0.1. 
The density might be even lower because, as discussed in §12.2.21 
the accretion of particles with an impact parameter higher than 
200 kpc will be considerable (and we are considering in our cal- 
culation only particles within b < 200 kpc). Lower densities are 
more in the range of possible values of intracluster density in 
small groups. However, the velocity of 10-20 km/s is somewhat 
low. For comparison, the velocity of the centre of the Galaxy 
with respect to the barycentre of the Local Group is 81 km/s 
[Braun & Burton 1999; eq. (3)]. Maybe in the first stages of the 
halo formation, when the galactic masses were much lower than 
MhaioCO = L2 X 10 12 M , the velocities of the galaxies with 
respect to the average IGM were also lower, since the gravita- 
tional potential responsible for such dynamics produced lower 
accelerations. A more developed study is necessary to clarify 
the question. 

In any case, statistically, it is possible to find a significant 
fraction of galaxies with low velocities with respect to the av- 



erage IGM. The present mechanism would apply to these few 
galaxies. 



Up to now, we have considered the standard model of galax- 
ies with massive dark matter halos. However, some authors have 
questioned the very existence of such halos since its evidence 
is weak on galactic scales (Battaner & Florido 2000; Sellwood 
& Kosowsky 2001; Evans 2001; etc.) and the rotation curves 
or satellite motions can be explained by alternative hypotheses 
such as magnetic fields (Battaner & Florido 2000) or MOND 
(Sanders & McGaugh 2002). If we consider a M mo (t) = 10 11 
M , necessary to explain the rotation curve of the Galaxy up 
to 15 Kpc (Honma & Sofue 1996) together with a disc + bar 
mass of 6 x 10 10 M , we would find from eq. d40l > that (we keep 
ffoo = 45°) 

F > 0.1 o v oo (km/s)e" h '~ (kra/s) pc O (kg/m 3 ) > 8.5 x 10~ 25 , (43) 

in which the necessary density is 6 times lower than in expres- 
sion d42l . 

If the halo mass were still lower than ~ 10 n M , we could 
reach the conclusion that most of the halo mass was accreted 
by this mechanism {F close to unity). Here, however, we have 
not included any alternative effects such as magnetic fields or 
MOND, which would be necessary to make compatible the ro- 
tation curves with this low value of the halo. 

If the mass of the halo were much higher (Conroy et al. 2005, 
give typical masses for the halo of the galaxies of ~ 5 x 10 12 
M Q ), the relative fraction F would be around 15% lower than 
in expression (l42l . although, of course, the amount of accreted 
matter is nearly proportional to the halo mass, so in absolute 
terms the number of accreted solar masses would be higher. 

2.6. Strength of the bar 

In the present paper, we have fixed the parameters of the bar as 
given in 82. 1.41 which are typical for a normal strong bar: Q m = 
5.4 x 10 4 (in units 10 6 M Q kpc 2 , as in Athanassoula 1992). A 
weak bar would have Q m ~ 10 4 (Athanassoula 1992). Since the 
energy loss is proportional to A and, A is proportional to the mass 
of the bar [eqs ©, ©], it is straightforward to deduce that within 
the linear regime of accretion (F « 1) F will be proportional 
to mass of the bar. Hence, a weak bar, typically with a mass five 
times lower than that used here, would reduce F by a factor of 
five. 

3. Other kinds of dynamical friction 

3.1. Spiral arms and warps 

As said in §|2] when the gravitational force varies with time, we 
get a non-conservative potential; that is, a potential in which the 
energy of the particle in this potential is not constant. Apart from 
the bar/triaxial bulge, there are other non-axisymmetric struc- 
tures in the galaxy, which, due to galactic rotation, produce this 
variation in the potential; for instance, the spiral arms or the 
warps in the disc. However, these structures normally have a 
mass much lower than that of the bar (the spiral arm of the Milky 
Way would have a ~ 10% of the disc mass according to the 
model of Bissantz & Gerhard 2002; and the mass of the external 
disc, where the warp is present [Lopez-Corredoira et al. 2002b] 
is also lower than 10% of the disc mass); therefore, their effect 



2.5.3. Mass of the halo 
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Fig. 9. Graphical representation of the hyperbolic trajectory of a 
particles which a approach a star (at the centre O). 



should be considerably lower, except perhaps in some cases with 
prominent spiral arms. These cases will not be explored in the 
present paper. 

3.2. Interaction with individual stars 

A way to produce changes in energy is in the interaction of the 
accreted particles with individual stars. Since individual stars 
move, the potential of the interaction is not steady, although 
it is only predominant with respect to the global potential of 
the whole galaxy when the particle approaches too close to 
the star. The effect is similar to the exchange of energy of a 
comet/spacecraft in an approach to a planet the Solar System. 

The velocity of a particle when it approaches a star of mass 
M, with respect to this star is 



Vdif = v - v„ 



(44) 



where v is the velocity of the particle with respect to the centre of 
the galaxy, and v* is the rotation(+peculiar) velocity of the star 
around the centre of the galaxy. The final velocity of the particle 
when it recedes from the star after the interaction is 



Vdif =v - V*, 



(45) 



where v' is the final velocity with respect to the centre of the 
galaxy. The star may be considered point-like and the path of the 
particle in the interaction with the star as a hyperbola with im- 
pact parameter b», which follows (see subsection 3.3 of Lopez- 
Corredoira et al. 2002a) 



tan — = 



y GM, 



(46) 



where y is the angle of deviation of the trajectory due to the 
interaction as represented in Fig. [9] 

The final velocity with respect to the star will be: 



dif 



= \Vdif\, 



v dif = Vdi f cos W e "-dif + v dif sin(y)e ft _ di/ , 



(47) 
(48) 



with the unit vector as specified in Fig. [9] and with respect to the 
centre of the galaxy: 



v =v' +v, 

dif 



(49) 



= (Vdif cos y + v. cos a)e v -dif + (v d if sin y + v* sin a)eb-dif, 

COS or = V t Vdif. 

Hence, 

v' = ^v 2 . f + v 2 + 2v dif v t cos(a - y), (50) 

while the initial velocity with respect to the centre of the galaxy 
is 



v = a/ V *'/ + v * + 2v dif v * cos(a). (51) 
This implies an increment of the energy of the particle 

AE - —m(v' 2 - v 2 ) = mvdifV f [cos(a - y) - cos(a')]. (52) 
2 

An approximation for small interactions (y << 1), using ( |46] i, 

(53) 



1 , ,7 7, „ • GM, 
AE = -m(v — v ) = 2mv* sin a 

2 b t v dif 



The particles trapped by the galaxy due to this interaction 
will be those with AE < -^mv^, (vco as the velocity of the par- 
ticle at infinite distance always is); that is, those particles with 
impact parameter 



K < 



4GM* sin |a|v» 



VdifV a 



(54) 



a < 



For M* — 1 Mq, v* = 200 km/s (the rotation as the main com- 
ponent), v M = 50 km/s, for an average case of a — -45° and 
Vdif ~ 100 km/s, the constraint is b* <~ 2 A.U. With the typical 
surface density in the galactic disc of 10 - 100 star/pc 2 (Lopez- 
Corredoira et al. 2002b), the probability of a particle having an 
impact parameter less than 2 AU is ~ 3 - 30 x 10~ 9 . This is just 
a rough estimate but the order of magnitude should be more or 
less correct. Therefore, the accretion ratio is much lower than 
the values of 7?2oo obtained in 32.21 and a much lower accretion 
(totally negligible) is expected to be produced in terms of the 
star/particle interaction. 

4. Comparison with observations: the effect of the 
bar on the halo mass 

It is not straightforward to check this effect on real galaxies. We 
would need to know the complete history of each galaxy, the 
epoch of formation of its bar and other characteristics which are 
not directly available. Nonetheless, we can analyse whether spi- 
ral galaxies with a prominent bar have different rotation-curve 
amplitudes from the same type of galaxies without bar. Indeed, 
this analysis cannot be a definitive proof of the present mecha- 
nism because there are other possible alternative explanations in 
terms of the interaction between halo and bar. For instance, halos 
might allow bars to become stronger (Athanassoula & Misiriotis 
2002; Athanassoula 2002, 2003). In any case, the present section 
might be a first step in the search of this effect. 
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In order to do this, we use first the data from the Mathewson 
et al. (1992) catalogue on rotation curves, which contains data 
for 1355 spiral galaxies. We subdivide the sample into 12 
groups: morphological type from +3 to +8 (equivalently from 
Sb to Sd) and each of them with or without bar (either SB/SAB 
or S/SA; according to the classification in this survey similar to 
the RC3 catalogue). This is a rough classification, because in- 
deed most of the spiral galaxies should have a bar to some de- 
gree (e.g., Knapen et al. 2000, Fathi 2004). From Mathewson 
et al. (1992), we take the maximum rotation velocity. This cata- 
logue includes data from optical H a rotation curves and HI pro- 
files. The results are given in Table Q] The table, although not 
conclusive, seems to show larger rotation velocities for earlier 
types (for +4, the average rotation velocity of barred galaxies is 
2.9cr higher than the non-barred galaxies) and vice verse for later 
types (for +7, the average rotation velocity of barred galaxies is 
3.1cr lower than the non-barred galaxies). The first thing might 
be explained in terms of the mechanisms proposed in this pa- 
per: throughout their lifetimes, bars help to accrete matter and 
this is more remarkable in earlier types because the bars are 
more prominent in these types. The hypothesis that halos allow 
bars to become stronger is also a good possible explanation; in- 
deed, the simulations in Athanassoula & Misiriotis (2002) and 
Athanassoula (2002, 2003) refer only to early-type disc galaxies, 
with no or hardly any gas, and this would explain why this effect 
is not necessarily seen in later types. The second fact is intrigu- 
ing: it might have to do with the history of the galaxy, in which 
these late-type galaxies without bars had one in the past, or with 
a major role of the spiral arms in non-barred galaxy. Also, er- 
rors in the classification of the galaxy type could be the cause of 
this different velocity in barred and non-barred galaxies (the bar 
could contribute to wrong visual classification). 

HI profiles with extension of the wings up to a 20% of the 
maximum trace the halo mass within larger radius, so the width 
of this HI profile (HW2QMhi in TableQ} is also considered. The 
results are similar to the results with the rotation curves for late- 
type spiral galaxies, and no trend is found for early-type ones. 

The same comparison between barred and non-barred galax- 
ies could be useful if made in galaxies in clusters. There 
are very few galaxies of Mathewson et al. (1992) classified 
as SB and "in cluster": according to the SIMBAD database 
(http://simbad.u-strasbg.fr/), only 5 of the 95 barred galaxies are 
known to be in clusters ('GiC'), so this catalogue is not useful for 
this purpose. We can use another catalogue with fewer galaxies 
in total (329 galaxies) but with a higher ratio of cluster members: 
Vogt et al. (2004). Results are shown in Table|2]for all the galax- 
ies, and in Table|2]only for galaxies in clusters. Nothing new is 
observed with respect to the Mathewson et al. (1992) analysis, 
which shows more or less the same trends. Analysis of barred 
galaxies in clusters (there are in total 23 for our analysis; not 
many, but more than the 5 galaxies in Mathewson et al. 1992) 
does not show a clearer trend of distinction with the non-barred 
galaxies: table[3] 

It would be better to obtain information on the total mass of 
halos rather than the mass within a few tens of kpc, as is the case 
for rotation curves, but this information is not available for a 
large sample of galaxies, and the errors in the estimates (through 
satellites orbit measures; e.g. Conroy et al. 2005) are inaccurate 
and not available for a large sample of individual galaxies to do 
statistics. 



5. Conclusions and discussion 

A new mechanism has been proposed in the interaction of bars 
with the surrounding medium, which produces a variation of en- 
ergy of the particles crossing the galaxy. On average, the energy 
gained is equal to the loss energy by these particles, so this does 
not produce variations of energy in the average halo. However, 
for the external IGM particles which cross the galaxy, this mech- 
anism causes some of them to be trapped to form part of the halo 
because their velocities are reduced to less than the escape ve- 
locity. By means of this, an extra amount of IGM matter can be 
accreted onto barred galaxies. The increase in the evaporation of 
the halo due to bar dynamical friction, however, is negligible. 

Rough calculations of the amount of accreted matter show 
that this must be negligible in most cases with normal expected 
conditions of density and galaxy-IGM velocity, and only in 
galaxies of galaxy-IGM velocities <x 70 km/s with IGM den- 
sity similar to that of the cluster of galaxies does it produce sig- 
nificant (> 10%) amounts of halo mass accreted by this mech- 
anism. The IGM density might be lower with motions of the 
galaxies close to the IGM average motion, or with motion per- 
pendicular to the plane of the galaxy. Clusters of galaxies also 
have other mechanisms of galaxy stripping that can counteract 
the accretion; indeed, the IGM and halos are mixed and we could 
even talk about a common halo for all the galaxies in a cluster. 
This means that we might find that there is no important increase 
in mass in the barred galaxies of clusters. 

The analysis of observational data of rotation curves shows 
no clear trends. There is a slight increase in the average mass 
of early-type barred spiral galaxies with respect to non-barred 
galaxies of the same type, and viceverse for late-type spiral 
galaxies, but some sort of selection factor might explain it. No 
effect could be observed for galaxies embedded in clusters, al- 
though we cannot exclude the possibility that the effect proposed 
in this paper is present to some non-negligible degree because 
it might be lower than the error bars. In Tables Q] [2] and [3] the 
relative error bars of the velocities are higher than 10%, which 
implies a mass error higher than 20%, so possible slight average 
variations of less than ~ 20% in mass would not be detectable. 

The fashion in galactic models nowadays is to include very 
massive halos whose matter comes from the evolution of the ini- 
tial fluctuations of the large scale structure (e.g., Betancort-Rijo 
& Lopez-Corredoira 2002) and the accretion of dwarf galax- 
ies, the hierarchical scenario (e.g., White & Frenk 1991). The 
mechanism presented in this paper is part of the possible mech- 
anisms of accretion and it can also be applied to the satellite 
dwarf galaxies. Nevertheless, if the very massive dark matter 
halos did not exist and the rotation curves were explained by 
mechanisms different from dark matter (e.g., MOND, Sanders 
& McGaugh 2002; magnetic fields, Battaner & Florido 2000), 
the present mechanism could be much more important as a rel- 
ative contributor to the accretion of matter onto the galaxy, and 
it would not be detected by analysing rotation curves because 
these would not be related to the halo mass. 

What is noteworthy in this paper is that the present mecha- 
nism does not suppose any speculative scenario about the galax- 
ies themselves. The mechanism should be present "always". The 
only question is whether the amount of accreted matter in barred 
galaxies is significant enough to be detected or not, because this 
depends on the conditions of the IGM. 

Acknowledgments: Discussions with J. E. Beckman (IAC, 
Tenerife, Spain) have produced some ideas which were used as a 
basis for the development of the present paper. Thanks are given 
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Table 1. Comparison of v max [km/s] (maximum velocity of rotation curves) and HWIQMhi [km/s] (Half of the HI profile between 
points where the intensity falls to 20% of the highest channels in each half of the profile and corrected for relativistic effects; Note: 
average only with the galaxies which have this datum available) for spiral galaxies in Mathewson et al. (1992) catalogue in barred 
and non-barred cases as a function of the galaxy type. 



Galaxy type Non-barred: (N), OwX (HW20M HI ) Barred: (N), <v max ), (HW20M HI ) 



+3 (339), 174 ±3, 164 ±5 (23), 182 ± 15, 180 ± 17 

+4 (114), 155 ±5, 156 ±7 (10), 205 ± 17, 158 ± 20 

+5 (96), 127 ± 6, 127 ± 9 (18), 132 ± 13, 125 ± 13 

+6 (526), 150 + 2, 150 ±3 (31), 127 ± 8, 116 ± 6 

+7 (26), 112 ±9, 117 ± 13 (8), 77 ± 7, 81 ± 7 

+8 (55), 91 ± 4, 96 ± 4 (11), 78 ± 5, 85 ± 5 



Table 2. Comparison of OW3 [km/s] (measured velocity width in optical profiles with corrections for the shape of the rotation 
curve, cosmological stretch and converted to edge-on viewing measured at twice the scale length), RW2 [km/s] (measured 21 cm 
line profile velocity corrected and converted to edge on viewing; Note: average only with the galaxies which have this datum 
available), M/L [solar units] (corrected mass to light [I-band] ratio; Note: average only with the galaxies which have this datum 
available) for spiral galaxies in Vogt et al. (2004) catalogue in barred and non-barred cases as a function of the galaxy type. There 
are no barred galaxies with available data in this catalogue of types +6 or later. 



Galaxy type Non-barred: (N), (OW3), (RW2), (M/L) Barred: (N), (OW3), (RW2), (M/L) 

+3 (94), 353 ± 12, 388 ± 14, 0.47 ± 0.03 (25), 404 ± 26, 434 ± 35, 0.44 ± 0.06 

+4 (37), 321 ± 16, 359 ±21, 0.52 ±0.04 (3), 419 ± 36, 391 ± 57, 0.76 ± 0.09 

+5 (83), 303 ±11, 331 ±10, 0.46 ±0.03 (9), 252 ± 28, 295 ± 43, 0.52 ± 0.11 



Table 3. Same as Table |2]but only with galaxies considered bona fide members of a cluster of galaxies. There is only one barred 
galaxy of type +4 "in a cluster" of Vogt et al. (2004) catalogue, which is not included. 



Galaxy type Non-barred: (N), (OW3), (RW2), (M/L) Barred: (N), (OW3), (RW2), (M/L) 
+3 (73), 350 ± 12, 393 ±15, 0.47 ±0.03 (19), 384 ± 31, 383 ± 32, 0.40 ± 0.07 

+5 (33), 317 ± 17, 346 ± 16, 0.49 ±0.05 (4), 234 ± 55, 321 ± 91, 0.46 ± 0.19 



to E. Athanassoula (Obs. Marseille, France) and the anonymous 
referee for helpful comments on the draft of this paper and sug- 
gestions to improve it; and T. J. Mahoney (IAC, Tenerife, Spain) 
for proof-reading of this paper. This research has made use of 
the SIMBAD database, operated at CDS, Strasbourg, France. 
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